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^ ; Abstract 

. The path-integral of the fermionic oscillator with a time-dependent frequency is ana- 

lyzed. We give the exact relation between the boundary condition to define the domain in 
' which the path-integral is performed and the transition amplitude that the path-integral 

calculates. According to this relation, the amplitude suppressed by a zero mode does not 
indicate any special dynamics, unlike the analogous situation in field theories. It simply 
says the path-integral picks up a combination of the amplitudes that vanishes. The zero 
mode that is often neglected in the reason of not being normalizable is necessary to obtain 
the correct answer for the propagator and to avoid an anomaly on the fermion number. 
(D [ We give a method to obtain the fermionic determinant by the determinant of a simple 

i-C 2x2 matrix, which enables us to calculate it for a variety of boundary conditions. 

> : 

PACS: 11.30.Fs; 03.65.-w; 02.70.Hm; 11.30.Pb 

.«?; 

Fermionic determinant of the operator D is the one we first encounter in the analysis of 
the quantum physics in the path-integral formalism. It is the exact result of the Grassmann 
path-integral made from a bilinear Lagrangian ipVip over the fermionic degrees of freedom if) 
and ip. The determinant carries an important information about the time evolution of the 
fermions under the influence of bosonic background. Especially when T> or its adjoint X>t 
has a zero mode, the zero-frequency eigenmode of the operator, the determinant vanishes and 
corresponding transition is suppressed. A typical example of such situation happens with D the 
Dirac operator in SU(2) gauge theory. It possesses zero modes in the instanton background. The 
consequent suppression of the transition is interpreted to reflect the fermion number violation 
due to the anomaly on the fermion current []]] . 
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The physical significance of the zero modes of Dirac operators was first advocated in Ref. [0 



and has been discussed in various places in physics. We still notice there seems to be a confusion 
on the treatment of zero modes in the case of fermionic oscillator, the simplest system containing 
only one fermionic degree of freedom. Its Lagrangian ipT>if) in imaginary time formalism is given 
by a simple first order differential operator 



with respect to the imaginary time r, where v(r) is the time-dependent angular frequency 
induced by a coupling to a bosonic degree of freedom. We also define 



Assume v(t) changes its sign along the evolution in the time interval [-T/2, T/2], for example, 
as v(t) = tanh(r). The solution of Dip = is then regarded as the zero mode since it is 
normalizable in the limit of T — > oo. Gildener and Patrascioiu have argued by an explicit 
calculation that there is no zero mode available in the determinant calculus even in this simple 
example ||, while Salomonson and Van Holten have taken advantage of the zero mode in their 
calculation for the supersymmetry breaking |4[]. 

In addition to the question of the existence of zero mode, we are concerned that the fermionic 
oscillator would have an anomaly if the numbers of the zero modes of D and are different. In 
the simple example mentioned above, the solution of T)'</> = is not thought to be the relevant 
zero mode: for sufficiently large T, it becomes zero almost everywhere when normalized to 
one. Thus there appears to be the asymmetry in the numbers. According to the path-integral 
formulation of the anomaly [||], this asymmetry induces a phase in the path-integral measure 
under the global phase transformation, ip — > e ip and ip — *> e~ ie ip, which indicates the non- 
conservation of the fermion number. This contradicts our naive intuition, that the fermionic 
oscillator should have no anomaly since we can calculate any amplitude with no regularization. 
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Motivated by this question of zero mode, we will clarify in this letter the exact calculation of 
the path-integral of the fermionic oscillator. In doing the path-integral, we first prepare an or- 
thonormal complete set for all possible configurations. The eigenmodes of a certain self-adjoint 
operator are known to constitute such a complete set. The useful choice for the calculation is 
to use the differential expressions V^V and as the candidate for the operator. As is known 
in mathematics, we need to impose the boundary condition to make the differential expressions 
to be self-adjoint. The important observation found in this letter is that the path-integral 
represents a different transition amplitude if the boundary condition is different. We will give 
the exact relation between the boundary condition and the transition amplitude, and verify it 
by comparing the path-integral results with those obtained by the operator formalism. The 
suppression for the zero mode will be understood within this boundary condition dependence. 

We first carry out the calculation in the operator formalism. The Hamiltonian of the 
oscillator is 

H(t) = v{t)\ (V^ - tftft) , (3) 

where \I/ and ^ are the annihilation and creation operators in the two-dimensional space 
spanned by vacant |0) and occupied |1) states, 

*|0> = * t |l)=0, tf|l) = |0), tf f |0) = (4) 



From this Hamiltonian, the evolution operator from the initial time Tj to the final Tf is obtained 
by 



U(T f ,T t )=Tex P 



dT'H(r') 



(5) 



where T represents the time-ordered product. The Hamiltonian commutes with the fermion 
number operator \J/t\l/, and it is obvious that the matrix elements of U are written as 



(0|C/(T / ,T,)|0)=exp 



1 fTf 



2J Tl 



drv(r) 



, (l|?7(7>,r i )|l>=exp 



l r T f 



2J Tl 



drvir) 



(6) 



and two off-diagonal ones equal to zero. These are all we have to know to obtain any transition 
amplitude. 



We turn to the path- integral formalism and consider 

cT/2 



[dip dip] exp 



dripVip 

-T/2 



(7) 



where the functional measure is defined by discretizing imaginary-time ||. We will come back 
this point later. Eq. (0) is related to the matrix elements of U(T/2, —T/2). Before seeing it, 
let us proceed the calculation of Eq. (|7|) usually adopted in literatures. The most decent way 
of the calculation for non-Hermitian T> such as the one in the present case is to choose the 
domains D v and of the square-integrable functions in the interval [—T/2, T/2], in which 
J)'T> and Z>P' are self-adjoint and positive semi-definite, and the non-zero eigenmodes have 
the one-to-one correspondence (see for example 0). The normalized eigenmodes (p( n \& D v ) 
and 0( n >(e D+) (n = 1,2,3,....), 

ptp^W = \ n ^ n) , £>£>V (n) = \ n (p {n \ (8) 

constitute a complete orthonormal set and they are related by 

1 -VtpW = _^x>V (n) = v {n \ (9) 



except zero modes. The eigenvalues so obtained are non negative. Using the expansions with 
these eigenmodes, 

iP(r) = £a n ^(r), ^(r) = £ «n0 (n) *(r), (10) 

n n 

we employ the integration measure [ete n rfa n ] instead of [dip dip] and obtain 

1/2 



I = M j [dada]e~ £ n ^ s " a " = TV [det (£>+£>) 
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where A/" is the Jacobian between the measures, det(X>'X>) is the infinite product of the eigenval- 
ues. The latter is well-defined in the combination with M . A different background v{r) results 
in a different complete set and different amplitudes in Eq. fllQ|), say b n and b n . The Jacobian 
between the measures [dacfa] and [dbdb] is the determinant of the matrix that represents the 
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linear transformation between the eigenmodes in the different sets. Since each set is complete 
and orthonormal anyway, the matrix is unitary and the Jacobian is a phase. Thus background 
dependence of M in (|TT|) can be at most a phase factor^. In the explicit examples adopted in 
later calculations, however, we can choose the eigenmodes as real function of r. This is because 
either the differential equations (H) or boundary conditions (see Eqs. (|26|) and (^7])) to solve the 
eigenvalue problem does not have any imaginary variable or constant. Thus the phase factor 
is 1 or —1. The Jacobian J\f then becomes background-independent except a possibility that it 
changes the sign at some special background. 

Although the differential expression V^V or VV^ seems self-adjoint by its form, it is not 
immediately so. One needs to define the domains and properly in the functional space 
by imposing the boundary condition. We notice that any ipi and y?2 in D v should obey 

^(-T/2)£Vi(-7V2) - Lp* 2 (T /2)Dlpi(T /2) = 0. (12) 

Under the usual definition of the inner product, 

fT/2 

(<f2,(fi)= drip* 2 (r)ipi(T), (13) 

J-T/2 

Eq. (0) means (tp2, T> j! T><pi) = (Uip2,T)<pi) which guarantees that T>^T> is non-negative, and 
leads to (ip 2 , V^Vipi) = (V^V(p2 1 y?i) which holds if it is self-adjoint. Note also that the equation 

fT/2 _ ,T/2 

/ dr^V^ = / drCD^M. (14) 

J-T/2 J-T/2 

holds as far as we expand ip and i/j as in fllOD with <p> and <fi that obey @. This equation means 
that any <p> G D v and any <fi G satisfy 

<f>*(-T/2)ip(-T/2) - 0*(T/2V(T/2) = 0. (15) 

Since Dtp G if ip G D v , Eq. ( |I~2D and (|15|) are equivalent. 

lr The procedure described here has a greater advantage in this property than the direct use of the eigenmodes 
of non self-adjoint T> as was done in Ref. @. 



In fact the mathematical theory of the self-adjoint extension of the differential expression 
tells us that we need to specify two linearly independent boundary conditions on the values, 
(p(—T/2) and <p(T/2), and the first derivatives, <p(—T/2) and <fi(T/2) in the present case of 
the second order differential expression. We refer the readers to the mathematical textbook J7| 
about the details. The problem is that the boundary condition so obtained is not unique. Then 
there naturally occurs a question about the boundary condition dependence of the path- integral. 

To see the boundary condition dependence, we review the derivation of Eq. (0) based on 
Ref. ||. We prepare the states 



|0) = |o) + |i)0, { e\ 



9{1 1 



(16) 



making use of Grassmann numbers 9 and 9. They satisfy the completeness relation, 



d9d9(l-99)\9)(9\ = |0> <0| + |1)(1| 



(17) 



The path-integral ([j]) calculates the evolution from \9 ) to \9 N ), (9n\U\9 ), where we have 
abbreviated U(T/2, —T/2) to U. To evaluate this, we discretize the time interval into N 
segments, each of which has the length e = T/N and write U as 



/ N 

U= lim T n(l-e%)) 

N— >oo \ , 
\n=l 
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where r n = nt — T/2. We insert (N — 1) pairs of Grassmann integrals of (9 n , 9 n )(n = 1, iV — 
1) in the form of the completeness relation flT^ ) as the (N — 1) junctions of the N factors 
(1 — eH{r ri )){n = 1, N). Using the explicit form (Q) and (|j), we obtain 



7JV 



U\0 ) = lim / de N - 1 dO N _ 1 ...d9 1 de 1 



' , v(r N ) ( v(t n ) \ s 

1 + Hi- e— — 9 N 9 N _ ± 



1 + - (l + £ %^) B N -^ + (l - J,.^' 



X... 



'i - ^ Mo 



(19) 
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All the bracket factors except the first one have the form 



, (i=l,...,N-l) 



(20) 



and can be safely replaced with 

exp <^ -e0j 



(21) 



because 9\ = and / d9i = 0. Then we arrive at a very close expression to Eq. (|7|) when we 
regard 9 n as ip{r n ) and 9 n as ip(r n ). Note, however, a difference of the first bracket factor in 
Eq. ([19]) from the others. We also wonder where 9n and 8 have gone when we got the result 

These questions are solved naturally by assuming that the exact definition of the measure 
[difidif)] in Eq. (0), corresponding to the way we have integrated it out, does include the integral 
over Qjy and 9 as 

I = J d6 N d6 N (l - 6 n 6 n )(6n\U\6 } (22) 

where 8q and On are related by the boundary condition to make Z>'Z> and T>T>^ self-adjoint. 
The relevant terms in the integrand in (p2|) is those quadratic in the Grassmann variables, 



;i - o N 9 N ){9 N \u\9 ) = ... - e N e N {o\u\o) + e N e {i\u\i) + ... 



This shows 



I=(0|[/|0)+/3<1|*7|1> 



(23) 



(24) 



if the boundary conditions imposed to define D v is 9 + 



'N 



0, or equivalently 



(p(-T/2) + /3(p{T /2) = 0. 



(25) 



Eqs. ( p4]) and (|25| ) are consistent with the fact that the fermionic path-integral presents the 
trace of U when carried out in the anti-periodic configurations {(3=1). 
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We verify Eq. ( P^j ) by explicitly calculating the determinant det(X^T>) at various boundary 
conditions. The boundary conditions we are interested in is written as 

<p(-T/2) + /MT/2) = 0, PV<P(~T/2) + V<p(T/2) = 0, (26) 

and 

P4>(-T/2) + 4>(T/2) = 0, X>V(-T/2) + /3X>V(T/2) = 0, (27) 

with a real parameter /3. This condition includes the periodic one at f3 = — 1 as well as the 
anti-periodic at 1. It is also available for the domain in which the zero mode lives. Once we 
put the boundary condition on (p, Eq. ([12]) and the one-one correspondence between f and </> 
determine the conditions on T>(p, <f), and T)*<j). The operators £>'X> and T?D* are proved to be 
self-adjoint under these boundary conditions JJJ. To calculate detiV^V), the 2x2 matrix 

M( A = ( tt!(z; -T/2) + p Ul (z; T/2) « 2 (z; -T/2) + /3u 2 (z; T/2) \ 
W " ^ -T/2) + T/2) /3D« 2 (2; -T/2) + D« 2 ( 2 ; T/2) J 

plays the central role, where Ui(z;r) (i = 1,2) are the linearly independent solutions of the 
equation 

V ] Vui(z-T) = zui(z;r) (29) 
and the parameter z is complex in general. We fix the normalization of these solutions by 

ui(*; -T/2) = 1 u 2 (z;-T/2) = 

u x {z\ -T/2) = ii 2 (z; -T/2) = 1. (30) 

Note the differential equation (|29|) does not have the first derivative term. Thus the Wronskian 
{u\Ui — Uiu 2 ) conserves and is equal to 1 at any r. The zero points of det M coincide with the 
eigenvalues of T>^T>: if A is one of the eigenvalues, there exists a non-trivial linear combination 
7iMi(A; r) + 72^2(^5 r) that satisfies Eq. (PB|); the equation for 7$ turns out to be 
M(X) i j'~fj = = 1,2) and thus detM(z) is zero at the eigenvalue A; one can reverse this 
argument in the opposite direction. 
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We then see the ratio [det(X^Z> — z)j det M(z)] is independent from the background v{r). 
The proof is essentially the same as the one given in Ref. || in the calculation of the other type 
of determinants. Let us consider two different operators T>\T>i and T> 2 T> 2 containing different 
backgrounds t>i(r) and v 2 (r), and denote their n-th eigenvalue by Ai ;n and A 2;n , respectively. 
Correspondingly let M\ and M 2 denote the matrix made by ( p8f ) and (^) with T>\T>\ and T> 2 T> 2 , 
respectively. The ratio defined by 

dgt(Pfoi - z) ^ ™ ( X 1;n -z \ 

det(p|p 2 -z) iivv-v 

is a meromorphic function of z, and it has a simple zero at each Ai ;n and a simple pole at each 
\ 2 - n . It goes to one as z goes to infinity in any direction except along the real positive axis. The 
ratio [det Mi(z)/ det M 2 (z)] is also a meromorphic function that has poles and zeros at exactly 
the same z. Note further we obtain 

det M(z) = Vu 2 (z; T/2) + (3 2 [ Ul (z; T/2) - v(-T/2)u 2 (z] T/2)} + 2(3 (32) 



using the condition (|30| ) and the conservation of the Wronskian. For sufficiently large 
^J\z\ 3> \v{t) 2 — v(t)\T, the frequency v(r) in (|29D becomes negligible. The solutions Ui(z;r) 
is then well- approximated by their free solutions, and ii 2 (T/2) ~ ui(T/2) ~ e^~* T /2. The 
determinant of M then grows as (1 + f3 2 )e^~* T /2 at sufficiently large \z\ (except along the real 
positive axis) independently from the backgrounds. The ratio [det M\ (z)/ detM2(z)] also goes 
to one in the same limit. Thus 

detjVlV! - z) = detMx(z) 
det(vlv 2 -z) det M 2 (z)' 

Eq. (|33|) establishes that [det^^V — z)j det M(z)] is background independent. 
We can now write?] 

M[det(V ] V)] l/2 = Af' [det M(0)] 1/2 . (34) 

2 Similar formulae that relate the determinant of differential operators with that of a matrix have been found 
in condensed matter physics H. 
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The factor Af' is background independent as is Af. The calculation of det M(0) is elementary. 
We obtain 

u 1 (0',t) = x 1 (t)+v(-T/2)x 2 (t), u 2 (0-t)=x 2 (t), (35) 

where 



exp 



dr'v(r') 



T/2 



yx{T) = exp 



dT'v(T') 



T/2 



/\\2 



x 2 (r)=x 1 (r) / dr'Mr')) 

-T/2 

(36) 



(37) 



(a 



Putting these solutions into (j32|) , we find 

detM(O) = ([l/i(r/2)] 1 / 2 + /3[xi(T/2)] 1/2 ) 2 . 

Note |/i = {xi)~ l . Recalling Eqs. (0), O), and (137|) , we obtain 

/ = AT'{[y 1 (T/2)] 1 / 2 + /3[x 1 (T/2)] 1 / 2 }. 

This is exactly Eq. (p^D, where the matrix elements are calculated explicitly by Eqs. (|6]) and 
(|3T)|). The factor A/ 7 turns out to be 1. We like to make a comment here. / in Eq. Q3"8] ) becomes 
negative, for example, when (3 = — 1 and x\(T/2) > yi(T/2). It means that we have chosen 
the negative solution in taking the square-root of det M in this case. This choice is justified by 
considering that I should be analytic with respect to the functional variation of v (r). The flow 
of V\ n that appears in (||) and (JTT|) must be smooth when the background changes continuously. 

We are prepared to answer the question of zero mode. Note that the zero mode candidate is 
X\ or yi in (|36|) . Let us start with assuming that the normalizable one in the usual sense is x\. 
The domain in which X\ resides is given by the boundary condition ( p6|) with (3 = —y\(T/2). 
Now that zero is an eigenvalue, the path-integral is zero. This is readily verified by putting the 



value of (3 into (38). The reason of the vanishing amplitude is, however, quite different from 
our interpretation of the same occurrence in the SU(2) gauge theory. There is no dynamical 
reason; we have just chosen a vanishing combination accidentally by the boundary conditions. 

We also notice the other zero solution y\ cannot be neglected. It satisfies the boundary 
condition Eq. ( p7|) with the same value of j3. Since we can normalize it any way as long as 
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the time interval T is finite (no matter how long it is), we cannot find any legitimate reason 
to abandon it. We can confirm the necessity of y\ in the calculation of the propagator. Let us 
define 

F(r,r') = 0(t -t')U(T/2,t)MJ(t,t')&U(t',-T/2) 

-6{t' - r)U(T/2, t')¥U{t', r) W(t, -T/2), (39) 

and consider its path-integral representation 

/•T/2 

^(r)^(r'). (40) 





fT/2 


j [dipdijj] exp 


— / dripVip 




J-T/2 



The exact relation of G integrated in the domain defined by Eqs. ( p6|) and (|27|) to the corre- 
sponding matrix element of F is obtained by applying the same argument that leads to Eq. (|21|). 
It is 

G{t,t') = 6(t-t') (0\U(T/2,t)W(t, t')&U(t', -T/2) |0) 

+0(r'-r)i/ 1 (T/2)(l|f/(T/2,r / )^ t t/(r / ,r)W(r,-T/2)|l). (41) 

In the operator formalism, Eqs. (|J) and (|) yield 

(0\U(T/2,t)W(t,t')¥U(t',-T/2)\0) = y 1 (T/2) 1 / 2 x 1 (T)y 1 (r'), 
(1\U(T/2,t)¥U(t',t)W(t,-T/2)\1) = x 1 (T/2) 1 / 2 x 1 (r)y 1 (r / ), (42) 

and, thus, using Eq. (|4I| ) we obtain 

G(r,T')=y 1 (T/2) 1 / 2 x 1 (T)y 1 (r'). (43) 



Interestingly the final result does not have any remnant of the time-ordered procedure in (|4l|) . 
We do the corresponding path- integral (0). The normalized zero modes are 

Mr) = T ;* (T) , ^>(r) = J' (T) . (44) 

H-rr- dT ' XliT "> 2 iL,^'^ 2 
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The Grassmann variables a x and di in the expansion (|H]), the coefficients of ip^ and <jP~\ do 
not appear in the action. Only the tp and ip in the integrand in ( [4*01) can supply them, and one 
gets 

G(t,t') =Af[det'(V^V)] 1/2 ip (1 \t)<J)V(t'), (45) 

where det'(P^P) is the product of the eigenvalues except zero. It is evaluated by 

det M(z) 



N 2 det'(V^V) = lim 



(46) 



For the calculation of the right-hand side in Eq. ("4T"), we use the fact that the solutions 
Ui(z; r){i = 1,2) have an expansion around z = 0, 



Ui(z; t) = Ui(0; t) + z 5ui(r) + ... 



(47) 



and the first order term is given by 



dUi(r) 



T/2 



dr' [ Ul (0; t)u 2 {0; t') - u 2 (0; r) Ul {0; r')] ^(0; r') 



(48) 



Since 5ui(—T/2) = 8ui(—T/2) = 0, the expansion pTj) satisfies the initial conditions ("30|). 
Putting Eqs. pTj) and ( flB]) into Eq. (|3~|) and using Eqs. (^) and (|36"|), we obtain 



T/2 
-T/2 



dry^rf 



fT/2 

/ dr xAt) 

J -T/2 



(49) 



Eqs. (["£"§), (P^D> (E§) and ( EH ) Si ye the correct result (fJ3|). We would have a wrong answer 
without </>W. The numbers of the zero mode belonging to D v and are the same independently 
of the specific boundary condition to define them. This is consistent with the absence of the 
anomaly in the fermionic oscillator. 

In summary, we have revealed and confirmed that the path-integral in different boundary 
conditions calculates different matrix elements of the time evolution operator. The exact rela- 
tion between the boundary condition and the matrix element tells us the reason of vanishing 
path-integral by the zero mode. It does not have a dynamical reason such as the fermion 
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number violation, but it occurs just because the zero mode takes the vanishing combination of 
the matrix elements. One may wonder whether the boundary condition dependence vanishes 
in the limit of T — > oo. It does not vanish. For example, the path-integral with the periodic 
boundary condition can be either positive or negative depending on the background, while that 
with the anti-periodic one is always positive no matter what the background is. It is natural to 
expect the existence of the similar boundary condition dependence in field theories. Then any 
conclusion based on a path-integral calculation should be spelled out with the boundary condi- 
tion, with respect to the imaginary time, explicitly specified. The conclusion might change for 
different boundary conditions. As far as we know there is little discussion on this dependence. 
A careful investigation of field theories taking consciously boundary conditions into account 
may reveal their new aspects. 
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